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Abstract 

Bats may be natural reservoirs for a large variety of emerging viruses, including mammalian coronaviruses (CoV). The 
recent emergence of severe acute respiratory syndrome-associated coronavirus (SARS-CoV) and Middle East respiratory 
syndrome coronavirus (MERS-CoV) in humans, with evidence that these viruses may have their ancestry in bats, highlights 
the importance of virus surveillance in bat populations. Here, we report the identification and molecular characterization 
of a bat p-Coronavirus, detected during a viral survey carried out on different bat species in the island of Sardinia (Italy). 
Cutaneous, oral swabs, and faecal samples were collected from 46 bats, belonging to 15 different species, and tested for 
viral presence. Coronavirus RNA was detected in faecal samples from three different species: the greater horseshoe bat 
(Rhinolophus ferrumequinum), the brown long-eared bat (. Plecotus auritus ), and the European free-tailed bat ( Tadarida 
teniotis). Phylogenetic analyses based on RNA-dependent RNA polymerase (RdRp) sequences assigned the detected CoV 
to clade 2b within betacoronaviruses, clustering with SARS-like bat CoVs previously reported. These findings point to the 
need for continued surveillance of bat CoV circulating in Sardinian bats, and extend the current knowledge on CoV ecology 
with novel sequences detected in bat species not previously described as p-Coronavirus hosts. 
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Introduction 

Bats (order Chiroptera) are one of the most diverse and 
widely distributed group of mammals, comprising more 
than 1300 species across 21 families [1], and play a key 
role in terrestrial ecosystems [2]. In the last few decades, 
they have been associated with several emergent zoonotic 
agents and recognized as potential reservoir of zoonotic 
viruses belonging to different families, such as Hendra, 
Nipah, Ebola, Lyssa, and SARS Coronavirus [3]. Bats rep¬ 
resent ideal hosts and reservoirs of viruses, due to their 
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high diversity and ecological characteristics that potentially 
facilitate virus transmission, such as the roosting together 
behavior and their ability to fly over distance [4, 5]. There¬ 
fore, it is essential to identify and characterize the currently 
unrecognized viruses circulating in bat populations. 

Coronaviruses (CoV, family Coronaviridae ) are envel¬ 
oped viruses with large single-stranded RNA genomes rang¬ 
ing between 27 and 32 kb [6, 7]. Based on genome-scale 
phylogenies, they are classified into four genera: Alphacoro- 
navirus (a-CoV) and Betacoronavirus (p-CoV), pathogenic 
for mammals, Gammacoronavirus (y-CoV) and Deltacoro- 
navirus (5-CoV), primarily (but not exclusively) circulating 
in birds [8, 9]. Most CoVs are disease-causing agents, asso¬ 
ciated with a range of respiratory, gastrointestinal, hepatic 
and neurological diseases both in humans and animals. They 
are also potentially transmitted among species, causing epi¬ 
demics after host transition [10-12]. Due to mutation and 
recombination events, new strains constantly originate and 
spread [4, 13]. 

Different bat species have been identified as natural reser¬ 
voirs of most mammalian a-CoV and p-CoV, including the 
Severe Acute Respiratory Syndrome SARS-CoV, and the 
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Middle East Respiratory Syndrome MERS-CoV [14, 15], 
which recently caused outbreaks in the human population 
[4]. SARS-related CoV was first identified in Rhinolophus 
bats from China [16], and subsequently many mammalian 
Co Vs were detected in Chinese bats [17, 18] and globally 
[19]. Further investigations identified a-CoV and P-CoV 
in several European countries from bats of various species 
[ 20 - 22 ]. 

In Italy, few studies have been conducted on bat CoV 
ecology: viral detection and prevalence was investigated in 
greater horseshoe bats [23], alpha- and betacoronaviruses 
were identified in different species, such as Kuhl’s pipist¬ 
relle, lesser horseshoe bats and serotine bats [10, 24], and 
recently two MERS related beta CoV detected in Italian 
Kuhl’s pipistrelle and Savi’s bat were fully sequenced [25]. 
A total of 21 bat species are currently recognized as liv¬ 
ing on the island of Sardinia (Italy), lying geographically 
at the interface between Europe and Africa, and recently 
indicated as a major Mediterranean hotspot for bat biodi¬ 
versity [26]. Nothing is known about CoV presence and 
distribution among bat populations on the island, and this 
is the first study aiming to fill this gap. Herein we report the 
molecular identification and typing of a SARS-like Corona- 
virus detected in three species of Sardinian bats belonging to 
Rhinolophus, Plecotus and Tadarida genera, during a viral 
survey conducted on non-invasive bat samples collected 
from different localities on the island. We also inferred the 
phylogenetic characteristics of the identified P-CoV, and its 
evolutionary relationships with other bat Co Vs identified in 
Europe and globally. 


Materials and methods 

Bat sample collection 

A total of 46 bats, belonging to 15 different species, were 
captured using mist nets or hand nets in 8 different locali¬ 
ties distributed in northern and central Sardinia: Grotta sa 
Rocca Ulari (Borutta), Galleria Casteldoria (Santa Maria 
Coghinas), Diga Santa Chiara (Ula-Tirso), Monte Coral- 
linu (Dorgali), Tilipera (Bonorva), Mularza Noa (Bolo- 
tana), Belvi, and Grotta Verde (Alghero), listed in Table 1. 
Sampling locations included caves, buildings, tunnels, 
woods, abandoned bridges, and other type of natural roosts. 
Captures were performed by expert and trained biologists 
within a regular census program of Sardinian bat popula¬ 
tions, between November 2015 and November 2016. Bats 
were placed individually into cotton bags before investiga¬ 
tion. Species, sex, age category, and forearm length were 
determined; oral and cutaneous swabs were collected using 
sterile micro-swabs (Aptaca), and then bats were released 
at their capture site. Faecal samples were readily collected 
from cotton bags or during bat handling and suspended in 
1 ml RNA later RNA stabilization solution (Qiagen) and 
then stored at - 20 °C before processing. 

DNA and RNA extraction, reverse transcription 

DNA was extracted from oral and skin swabs tips using 
the QiAmp DNA Mini Kit (Qiagen), following the 


Table 1 List of bat samples 
analyzed: host species, sampling 
locations, and type of samples 
collected 


Host species 

Sampling location 

Skin swab 

Oral swab 

Faeces 

Rhinolophus mehelyi 

1 

2 

— 

2 

Rhinolophus ferrumequinum 

1,2,8 

3 

3 

3 

Miniopterus schreibersii 

1 

9 

— 

7 

Myotis punicus 

1,2,5 

8 

8 

6 

Myotis capaccinii 

2,8 

3 

4 

3 

Myotis emarginatus 

2,6 

3 

2 

2 

Plecotus sardus 

3 

— 

1 

— 

Hypsugo savii 

4,6 

1 

1 

— 

Pipistrellus pipistrellus 

5 

2 

— 

— 

Pipistrellus kuhlii 

6 

— 

— 

1 

Plecotus auritus 

7 

1 

1 

1 

Plecotus austriacus 

8 

— 

1 

— 

Nyctalus leisleri 

6 

1 

1 

— 

Eptesicus serotinus 

8 

— 

— 

1 

Tadarida teniotis 

8 

1 

1 

1 

Tot 


34 

23 

27 
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Locations: (1) Grotta sa Rocca Ulari (Borutta), (2) Galleria Casteldoria (S.M. Coghinas), (3) Diga S.Chiara 
(Ula-Tirso), (4) Monte Corallinu (Dorgali), (5) Tilipera (Bonorva), (6) Mularza Noa (Bolotana), (7) Belvi, 
(8) Grotta Verde (Alghero) 
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manufacturer’s protocol specific for swabs. DNA was stored 
at - 20 °C for subsequent amplification analysis. After sus¬ 
pending faecal samples in phosphate buffered saline (PBS), 
RNA was extracted from faeces by using the QiAmp Viral 
RNA Mini Kit (Qiagen), according to the instructions sug¬ 
gested by the manufacturer. Extracted RNA was stored at 
- 80 °C. First-strand cDNA was synthesized utilizing the 
Affinity Script Multi Temperature cDNA Synthesis Kit (Agi¬ 
lent Technologies), following the random primers strategy 
protocol, and then stored at — 20 °C. 

PCR and RT-PCR, screening and sequencing 

DNA extracted from bat oral swabs was tested for the pres¬ 
ence of Herpesvirus using a universal nested PCR [27]. 
DNA extracted from skin swabs was screened for Poxvirus 
and Papillomavirus presence, by PCR based on pan-pox 
primers [28] and on primers targeting a conserved region 
ofPPVLl gene [29]. 

RT-PCR was performed to detect Coronavirus in bat 
faecal samples, by using two pairs of primers targeting a 
conserved fragment of RNA-dependent RNA polymerase 
gene (RdRp): 11FW and 13RV [23], and Corona 1 forward 
and Corona 2 reverse [10, 17, 30], targeting 170 bp and 
440 bp respectively. In a typical PCR, 2 pi of cDNA was 
amplified in a 50 pi reaction mixture containing 0.2 mM 
of each dNTPs, 0.2 pM each of Forward and Reverse prim¬ 
ers, IX TaqBuffer and GoTaq (Promega, Italy) DNA poly¬ 
merase. The cycling conditions were as follows: 5 min at 
95 °C followed by 40 cycles of 1 min at 95 °C, 1 min at 
50 °C, and 1 min at 72 °C, followed by a final extension 
time of 5 min at 72 °C. From the first amplified sequence 
of approximately 170 bp, in order to further amplify the 
RdRp gene, 4 sets of primers were newly designed, after 
aligning the sequence detected in this study with the most 
similar sequences obtained by BEAST analysis, and choos¬ 
ing conserved regions to design consensus walking prim¬ 
ers. The first 170 bp amplicon was extended by design¬ 
ing and using downstream primers within the fragment 


and upstream primers based on alignments of phyloge- 
netically related CoV from GenBank. Table 2 lists the 
primers designed in this study (BatCVIF/R, BatCV2F/R, 
BatCV3F/R, and BatCV4F/R). PCR products were purified 
using DNA Clean & Concentrator purification kit (Zymo 
Research, CA, USA), and directly sequenced in both for¬ 
ward and reverse directions by BMR Genomics Sequenc¬ 
ing Service (Padova, Italy). Forward and reverse sequences 
were edited and aligned with MUSCFE Multiple Sequence 
Alignment tool (http://www.ebi.ac.uk/Tools/msa/), align¬ 
ments were edited with BioEdit software [31]. Sequences 
were checked against the GenBank database using nucleo¬ 
tide BEAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi). 

Phylogenetic analysis 

Phylogenetic trees were estimated using the Maximum 
Fikelihood (ME) and Maximum Parsimony (MP) meth¬ 
ods implemented in MEGA 6 [32], after testing the best 
fitting evolutionary model to the data. Bootstrap support 
values and reliability of internal branches were calculated 
from 1000 replicate trees. Models with the lowest BIC 
scores (Bayesian Information Criterion) are considered to 
describe the substitution pattern the best. Non-uniformity 
of evolutionary rates among sites was modeled by using 
a discrete Gamma distribution (+G) with 5 rate catego¬ 
ries and by assuming that a certain fraction of sites are 
evolutionary invariable (+1). After generating a gap-free 
nucleotide alignment of 1090 bp containing the novel virus 
as well as of Coronavirus RdRp reference strains, the most 
likely nucleotide substitution model selected to analyze 
these data was TN93+G (Tamura-Nei, Gamma distrib¬ 
uted), while analyzing an alignment of a higher number 
of CoV RdRp sequences of 396 bp, the model T92+G was 
selected (Tamura 3-parameter, Gamma distributed). The 
CoV sequence obtained in this study was submitted to the 
GenBank database under accession number MG978754. 


Table 2 Primer sequences used 
and designed in this study for 
bat Coronavirus detection 


Name 

Target sequence 

Primer sequence 

Source 

11FW 

13RV 

RdRp 

5 '-TGATGATGCCGTCGTGTGCTACAA-3' 

5 '-TGTGAGCAAAATTCGTGAGGTCC-3' 

Balboni et al. [23] 

BatCVlF 

BatCVIR 

RdRp 

5 '-TACAGTGATGTAGAAACTCC-3' 
5'-AAGGTCTGTCTCAGTCCAGCA-3' 

This study 

BatCv2F 

BatCV2R 

RdRp 

5 '-TGCTGGACTGAGAC AGACCTT-3' 

5 '-GTGTTATC ATTAGTTAGC ATT-3' 

This study 

BatCV3F 

BatCV3R 

RdRp 

5'-CTAAGCGCAATGTCATCCCT-3' 

5'-GCATGTTAGGCATGGCTCT-3' 

This study 

BatCV4F 

BatCV4R 

RdRp 

5 '-GCATCAAGAACTTTAAAGCAG-3' 
5'-GAATACATGTCTAACATGTG-3' 

This study 
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Results 

Coronavirus RNA was detected in three out of 27 faecal 
samples: 11%, close to previously reported frequency [10], 
less than that observed in other studies [17, 30], unlike 
results obtained in central-southern Italy [33]. Positive 
samples were collected in different locations of northern 
and central Sardinia from three different bat species: the 
greater horseshoe bat (Rhinolophus ferrumequinum), the 
brown long-eared bat (. Plecotus auritus), and the European 
free-tailed bat ( Tadarida teniotis). Table 3 lists samples 
positive for CoV detected in this study, including the Gen- 
Bank accession number assigned to the novel sequence. No 
Herpesvirus, Poxvirus and Papillomavirus were detected 
in the oral and/or cutaneous samples collected within this 
Sardinian bats viral survey. 

After detecting CoV in faecal samples, a total of 
1110 bp of RdRp gene were successfully amplified and 
directly sequenced. The same betacoronavirus sequence 
was detected in all three positive bat species. Upon 
BLASTN analysis, the viral sequence identified was found 
highly similar to Bat coronavirus BM48-31/BGR/2008 
(GUI90215), detected in a Blasius’s horseshoe bat (Rhi- 
nolophus blasii) in Bulgaria [14], with 100% query cover¬ 
age and 96% identity. The obtained sequence also showed 
high nucleotide similarity (89%) with Bat SARS CoV 
Rs672/2006 (FJ588686) identified in a Chinese rufous 
horseshoe bat ( Rhinolophus sinicus) [34], and with Bat 
SARS-like coronavirus RsSHC014 (KC881005) from the 
same bat species [35]. A first phylogenetic analysis was 
performed based on 1090 bp amplified fragment of the 
RdRp gene including the CoV identified in Sardinian bats 
and other 16 sequences representatives of CoV genetic var¬ 
iability. The phylogram (TN93+G, Fig. 1) clearly showed 
that Sardinian bat CoV belong to the p-CoV genus, and 
clustered together with the Bulgarian bat CoV identified in 
a Rhinolophus blasii (GUI90215). To have a deeper look 
into clade structuring within genera and further character¬ 
ize the phylogenetic relationships of the detected CoV, 
especially with other European bat coronaviruses, another 
phylogenetic analysis was performed based on 396 bp of 


the RdRp gene, and including a total of 30 sequences. Fig¬ 
ure 2 shows the (T92+G) tree obtained, clearly assigning 
the Sardinian bat CoV to clade 2b within the p-CoV genus. 
The branching also showed that the detected sequence was 
closely related to a betacoronavirus previously identified 
in a greater horseshoe bat (Rhinolophus ferrumequinum) 
in Italy (BtCoV/Rhi_fer/Itl7/ITA/2009, GenBank number 
KC633199), as well as other bat CoVs detected in Rhi¬ 
nolophus spp. from Slovenia (GQ409747) [36], and again 
from Italy (i.e., KF500954) [10]. All other clade groupings 
were consistent with previous studies [14]. 

Discussion 

At least 21 bat species are currently recognized as living 
on the island of Sardinia, Italy [37]. Some species are com¬ 
mon and adaptable, inhabiting both wild and urban habitats, 
while others are rare and restricted to limited areas of the 
island. Recently an island endemic species was character¬ 
ized (Plecotus sardus) and an endemic lineage of Plecotus 
auritus was also identified [38]. The existence of endemic 
bat taxa on Sardinia highlights the island importance in the 
conservation of the European bat community. Some Sardin¬ 
ian bat populations appear genetically differentiated and iso¬ 
lated from any continental gene pool and should therefore be 
considered as an evolutionarily significant unit (ESU) [39]. 

This is the first report of Coronavirus detection in faeces 
from bat species in Sardinia. A SARS-like betacoronavirus 
was detected in three species (Table 3): the greater horse¬ 
shoe bat (Rhinolophus ferrumequinum), common cave spe¬ 
cies distributed all over the island; the brown long-eared 
bat (Plecotus auritus ), rare species inhabiting few locali¬ 
ties (especially mountain woods) in central Sardinia; and 
the European free-tailed bat (Tadarida teniotis ), common 
species usually colonizing rocky slopes and Mediterranean 
bushes as well as buildings and bridges in urban areas. 

The detected virus (SarBatCoVl) belonged to clade 2b, 
including most of the SARS-like CoVs found in bats, as 
well as human SARS Coronavirus (Figs. 1, 2). SarBatCoVl 
appeared related to another betacoronavirus previously iden¬ 
tified in the Italian greater horseshoe bat (BtCoV/Rhi_fer/ 


Table 3 List of bat positive samples for Coronavirus, and GenBank accession number 


ID 

Host species 

Sample type 

Sex 

Location 

CoV strain 

GenBank 

accession 

number 

15F 

Greater Horseshoe bat R. ferrumequinum 

Faeces 

F 

2 

SarBatCoVl 

MG978754 

3 OF 

Brown Long-eared bat P. auritus 

Faeces 

nd 

7 

SarBatCoVl 

— 

31F 

European Free-tailed bat T. teniotis 

Faeces 

F 

8 

SarBatCoVl 

— 


Locations (ref. Table 1): 2. Galleria Casteldoria (S.M. Coghinas), 7. Belvi, 8. Grotta Verde (Alghero) 
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60 

92 

51 


100 


r Bat SARS CoV Rs672/2006/R.sinicus/FJ588686 
- Bat SARS corona virus Rp3/2004/DQ071615 
r - SARS-related bat CoV/R.sinicus/KF294440 
Coronavirus Lyral 1/R.affinis/KF569996 
f- Bat SARS-like CoV RsSHC014/R.sinicus/KC881005 


99 


100 


50 


87 


69 


- BtRs-BetaCoV/YN2013/R.sinicus/KJ473816 
■ SARS-like CoV WIV16/R.sinicus/KT444582 

— Bat SARS coronavirus Rm1/R.macrotis/DQ412043 
Bat CoV isolate Jiyuan-84/R.ferrumequinum/KY770860 


100 


100 i- Bat SARS coronavirus Rf1/R.ferrumequinum/DQ412042 

— SarBatCoVI 

- Bat CoV BM48-31/BGR/2008/R.blasii/GU190215 
Bat coronavirus HKU9-3/EF065515 


Bat coronavirus HKU9-4/EF065516 


Bat coronavirus HKU4-1/EF065505 

-Bat coronavirus HKU6/DQ249224 


100 


Bat coronavirus BtCoV/512/2005/DQ648858 


] 


Beta 


Alpha 


0.1 

Fig. 1 Phylogenetic analysis based on 1090 bp of the RdRp gene. 
RNA-dependent RNA polymerase-based phylogeny including the 
novel bat coronavirus detected in this study (SarBatCoVI) and other 
16 RdRp sequences from bat coronaviruses. Coronavirus sequence 

Itl7/ITA/2009, GenBank number KC633199), and also 
clustered with a p-Coronavirus isolated from a Blasius’s 
horseshoe bat in Bulgaria (BM48-31/BGR/2008, GenBank 
number GUI90215). As shown in Fig. 2, the two clades 
within p-Coronaviruses are strongly supported in the phylo- 
gram, and clade 2c clustered many other bat Co Vs identified 
in Europe (Italy and Spain) from different species. These 
findings confirm a correlation between SARS-like CoVs and 
bats belonging to Rhinolophus spp., consistently with previ¬ 
ous studies [14, 40], although the same RdRp sequence, in 
this study, was also detected in two other species belong¬ 
ing to the genera Plecotus and Tadarida. As previously dis¬ 
cussed, phylogenetic analyses providing insights into viral 
transmission among bat species have revealed a general clus¬ 
tering by geographic location rather than by bat species [4] . 
However, our results also show that similar coronaviruses 
(with nucleotide identity up to 96%) are found in bats dis¬ 
tributed in geographically distant regions, such as the island 
of Sardinia, mainland Italy and Bulgaria. 

Since the identification of SARSr-Rh-BatCoV in bats 
of the Rhinolophus genus [16], an intense effort was dedi¬ 
cated to investigate CoVs in bats, and many SARS-like 
viruses were detected and characterized in various species 
worldwide [4, 19]. Recent studies also provided evidence 
that bats may be the original ancestor hosts for SARS and 


identified in this study is shown in bold and assigned to the P-CoV 
genus. Bootstrap percentage values from 1000 resamplings are 
located at nodes of the tree 


MERS viruses, which recently caused pandemic and heavy 
consequences in humans [41-43]. Therefore, it is para¬ 
mount to investigate and characterize CoVs distribution 
and diversity in bat species, before their potential emer¬ 
gence in humans. 

Full-length genome analysis of the detected CoV should 
be attempted to thoroughly understand phylogenetic rela¬ 
tionships between SarBatCoV 1 and betacoronaviruses previ¬ 
ously described in European bats. Further studies including 
larger number of samples as well as continued surveillance 
on viruses circulating in Sardinian species are also needed 
to understand the role played by bats in the maintenance 
and transmission of beta CoVs. Monitoring and research 
on potential zoonoses are essential, and the prevalence 
and distribution of CoVs in Sardinian bats should be fur¬ 
ther investigated in urban and rural environment. However, 
without specific evidence, bats should not be considered as 
a major concern for public health, especially given the key 
ecosystem services they perform [44]. Bats play essential 
ecological roles in mosquito and pest control, as pollinators, 
and are keystone species in natural ecosystems. On top of 
this, bat species are often protected by national and interna¬ 
tional laws, as their conservation status may be vulnerable 
or threatened, and Sardinia was recently indicated as a major 
Mediterranean hotspot for bat biodiversity [26] . 
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100 


98 


82 


97 


96 


KJ473819 BtMd-BetaCoV/HuB2013 
KY770860 BatCoV isolate Jiyuan-84 Rhi_fer 
— DQ022305 Bat SARS coronavirus HKU3-1 


85 


KT444582 Bat SARS-like CoV WIV16 
NC004718 SARS coronavirus 


97 L KC881005 Bat SARS-like CoV RsSHC014 


96 


89 

55 


80 


r SarBatCoVI 

L KC633199 Betacoronavirus BtCoV/Rhi_fer/lt17/ITA/2009 
r GU190215Bat coronavirus R.blasii BM48-31/BGR/2008 
L GUI 90230 Bat coronavirus Rhijer BM48-35/BGR/2008 
p GQ404797 Bat coronavirus SL01A0082/2008/SVN 
r KF500954 Betacoronavirus R.hipposideros/ltaly/243585/2012 
64 L KF500953 Betacoronavirus R.hipposideros/ltaly/196814/2011 J 
-DQ249214 Bat coronavirus HKU4 isolate HKU4-1 


2b 


99 

50 


DQ249217 Bat coronavirus HKU5 isolate HKU5-1 


96 


78 


V5 


70 


p KF500941 Betacoronavirus N.noctula/ltaly/206679-5/2010 
HQ184062 Betacoronavirus E.isa/M/Spain/2007 
- KC243390 Bat coronavirus BtCoV/8-724/Pip_pyg/ROU/2009 


JX869059 Human betacoronavirus 2c EMC/2012 


1 — HQ184059 Betacoronavirus H.sav/J/Spain/2007 


100 


46 


97 

55 


KF500946 Betacoronavirus P.kuhlii/ltaly/206645-53/2011 
KF500950 Betacoronavirus P.kuhlii/ltaly/330375-15/2012 
KF500940 Betacoronavirus H.sav/ltaly/206645-40/2011 
KF500944 Betacoronavirus P.kuhlii/ltaly/206645-29/2011 


2C 


Beta 


DQ249235 Bat coronavirus HKU2 isolate HKU2-1 

-DQ249224 Bat coronavirus HKU6 isolate HKU6-1 

HQ184058 AlphacoronavirusP.kuh/lprima/Spain/2007 
HQ184051 Alphacoronavirus N.las/C/Spain/2007 
KF500945Alphacoronavirus P.kuhlii/ltaly/206645-41/2010 
KF500945 Alphacoronavirus P.kuhlii/ltaly/206645-41/2010 J 


97 


60 


100 


Alpha 


0.2 


Fig. 2 Phylogenetic analysis based on 396 bp of the RdRp gene. Phy¬ 
logenetic tree of the partial RNA-dependent RNA polymerase gene 
(396 bp) of coronavirus strains found in bats and humans, showing 
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